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ABSTRACT  * 

It  is  shown  that  the  repulsive  (not  merely  the  positive) 
portion  of  the  Lennard-Jones  potential  quantitatively  dominates 
the  equilibrium  structure  of  the  Lennard-Jones  liquid.  A  simple 
and  accurate  approximation  for  the  radial  distribution  function  at 
high  densities  is  presented. 
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It  is  often  remarked  that  the  repulsive  forces  in  a  simple  liquid 

play  an  important  role  in  forming  the  equilibrium  structure  of  such  a 
12  3 

fluid.  ’  *  By  studying  the  spectrum  (Fourier  transform)  of  the  two-particle 


correlation  function, 


h(k)  =  p  J [g(r)-l]e  i~*E  dr 


(1) 


where  g(r)  is  the  usual  radial  distribution  function  and  p  is  the  average 
number  density,  we  are  able  to  formulate  this  idea  in  a  precise  and  useful 
way. 

In  order  that  a  convenient  basis  for  comparison  be  available,  the 

Lennard-Jones  fluid  is  considered:  a  classical  N-particle  system  in  which 

the  total  potential  energy  is  a  sum  of  the  Lennard-Jones  pair  potentials, 

w(r)  *  4e[(o/r)12  -  (o/r)^].  Our  attention  is  focused  on  h(k)  for  the 

Lennard-Jones  liquid.  In  particular,  for  liquid  densities  (po32.4)  we 

hypothesize  and  then  verify  the  following  statements:  1.  At  intermediate 

and  large  wave-vectors  Ckc^-rr) ,  the  quantitative  behavior  of  h(k)  is  dominated 

py  the  repulsive  forces  (the  attractive  forces  are  primarily  manifested  in 

the  small  wave-vector  portion  of  the  spectrum);  2.  For  high  densities 
3  A 

(po  £.65)  the  behavior  of  h(k)  even  at  small  wave-vectors  (kcg:r)  is  at 
least  qualitatively  determined  by  the  repulsive  forces.  Physically,  the 
first  statement  is  understood  once  it  is  recalled  that  h(k)  +  1  represents 
the  linear  response  of  the  fluid  structure  to  a  disturbance  of  wavelength 
2n/k.  While  a  short  wavelength  disturbance  will  probe  both  the  repulsive 
and  attractice  forces  in  a  fluid,  it  is  reasonable  that  the  presence  of  the 
harsh  repulsions,  rather  than  the  slowly  varying  longer  ranged  attractions, 
will  dominate  the  response  of  the  fluid  to  such  a  disturbance.  The  second 


statement  is  equivalent  to  asserting  that  the  correlations  in  a  simple 

liquid  are  almost  entirely  due  to  excluded  volume  effects  when  the  density 

is  high.  While  this  latter  part  of  our  hypothesis  has  been  stated  quali- 
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tatively  by  many  authors,  *  *  its  quantitative  validity  has  not  been 
appreciated. 

To  see  how  the  repulsive  portion  of  the  pair-potential  affects  the 
liquid  structure,  it  is  convenient  to  consider  a  reference  system  in  which  the 
pair-potential,  UQ(r),  contains  all  the  repulsive  forces  implied  by  the 
Lcnnard-Jones  potential,  w(r),  and  no  attractive  forces.  With  the  additional 
condition  that  the  reference  system  has  a  well  defined  thermodynamic  limit, 
uQ(r)  is  defined  uniquely: 

uQ(r)  *  w(r)  +  e,  r<2*^o 
.  r^o 

It  should  be  noted  that  the  repulsive  part  of  w(r)  is  the  portion  with  negative 

slope  and  not  the  portion  which  is  merely  positive.  In  terms  of  the  radial 

distribution  function  for  the  reference  system,  gQ(r),  and  its  spectrum, 

A 

h  (k) ,  the  hypothesis  stated  above  can  be  expressed  as  follows:  At  all 

A  A 

liquid  densities,  the  approximation  hQ(k)  =  h(k)  is  very  accurate  for 
ko^it;  at  high  densities  (po^.65),  g(r)  closely  resembles  gQ(r).  Thus,  the 

A 

task  of  verifying  our  hypothesis  is  reduced  to  determining  g0(r)  and  hQ(k). 
Rather  than  perform  expensive  machine  calculations  to  obtain  this  information, 
we  have  invented  an  approximation  method  to  describe  the  reference  system. 

Tests  have  been  performed^  which  indicate  that  our  treatment  is  certainly 
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accurate  enough  for  the  purposes  of  this  work.  To  discuss  the  method,  we 
introduce  the  function  y£)(r)  which  is  defined  by  the  equation 

g0(r)  =  yo(r)exp[-6uo(r)]  (3) 

where  6  is  the  reciprocal  of  Boltzmann's  constant  times  the  temperature. 
Physically,  yo(r)  gives  the  correlations  that  exist  in  the  reference  system 
beyond  the  range  of  the  reference  interaction,  uQ(r).  Since  this  interaction 
is  harshly  repulsive,  it  seems  probable  that  yQ(r)  can  be  approximated  by 
the  similar  function  appropriate  to  a  hard-sphere  system  of  diameter  d, 
yj(r).  For  this  reason,  we  consider  the  following  approximation: 

g0(r)  =  yd(r)exp[-BuQ(r)]  (4) 


or 


A  A 

h  (k)  =  h 
0  c 


00 


p/^  yd(r)^ 


-6u 


(4') 


where  h^Ck)  denotes  the  spectrum  for  the  hard-sphere  system  of  diameter  d, 

A 

and  ud  denotes  a  hard-core  repulsion  of  diameter  d.  By  recalling  h  (k)'s 
role  as  a  linear  response  function,  it  is  seen  that  Eq.  (4)  is  equivalent 
to  assuming  that  for  intermediate  and  long  wavelength  disturbances,  the 
response  of  the  reference  system  is  accurately  given  by  the  response  of  a 
hard-sphere  system.  Thus,  by  equating  the  long  wavelength  (small  k) 
responses  of  the  reference  and  hard-sphere  systems,  we  have  a  physically 
reasonable  density  and  temperature  dependent  criterion  for  determining  d: 


(5) 
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Finally,  for  lack  of  better  analytic  or  extensively  tabulated  hard-sphere 

functions,  we  use  the  analytic  solution  of  the  Percus-Yevick  equation  to 
a  6 

give  us  hj(k)  and  yj(r)  (o<r<2d) .  Though  it  is  fairly  accurate,  it 
represents  our  major  source  of  error  in  the  computations  reported  below. 

When  Eq.  5  is  solved  for  d(f3,p)  and  plotted  along  several  isochores, 
graphs  are  obtained5  that  are  very  similar  to  those  given  by  Vcrlet  (Fig.  17, 

Ref.  1)  for  the  effective  hard-sphere  diameter  needed  to  empirically  "fit" 

a  a 

a  hard-sphere  spectrum,  h^(k),  to  the  Lcnnard-Jo nes  spectrum,  h(k),  at 

intermediate  wave-vectors.  Thus,  we  believe  Vcrlet  was  observing  the  behavior 

we  have  hypothesized:  for  ka;>Ti,  the  spectrum  is  dominated  by  the  repulsive 

A 

(reference)  forces.  In  Fig.  1,  we  plot  hQ(k)  (as  obtained  from  our  treat¬ 
ment  of  the  reference  system)  at  two  thermodynamic  states,  (labeled  I  and  II) 
and  compare  with  Verlct's  molecular  dynamics  results.*  For  state  I,  the 
density  is  sufficiently  low  to  illustrate  the  dramatic  discrepancy  and  then 

A 

similarity  between  hQ(k)  and  the  true  I.cnnard-Joncs  spectrum  as  we  pass 
from  small  to  large  wave-vectors.  Incidcntly,  for  state  I,  the  value  of 
d  obtained  from  Eq.  5  is  (1.005)o;  Verlet  reports  a  hard-sphere  model  d  of 
(1.0)o.  From  the  curve  for  state  II,  we  see  that  at  high  densities,  the 
excluded  volume  effects  dominate  the  correlations  to  such  an  extent  that 
hQ(k)  is  also  in  good  agreement  with  the  Lcnnard-Joncs  spectrum  for  small 
wave-vectors.  The  value  of  d  for  this  thermodynamic  state  is  ( 1 . 025 ) u  ; 

Verlet's  hard-sphere  model  d  is  (1.026)o. 

A  A 

The  agreement  between  hQ(k)  and  h(k)  that  is  observed  at  high  densities 
implies  that  gQ(r)  should  closely  resemble  g(r)  (except  for  very  large  r)  at 
high  densities.  In  Fig.  2  a  plot  of  gQ(r)  as  described  in  fcq.  4  is  compared 


♦ 
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with  "exact"  Lennard-Jones  g(r);  the  thermodynamic  state  is  similar  to 

state  II  in  Fig.  1.  The  agreement  between  gQ(r)  and  g(r)  is  excellent.  To 

emphasize  how  accurate  it  is,  we  have  also  plotted  the  g(r)  obtained  by 

6a  7 

solving  the  Percus-Yevick  equation  for  the  Lennard-Jones  fluid.  *  We 

A 

mention  that  the  Fourier  inversion  of  hQ(k)  at  a  moderate  density  (e.g. 
the  thermodynamic  state  I  considered  in  Fig.  1)  agrees  poorly  with  the 
Lennard-Jones  g(r)  for  all  r .  This  behavior  is  a  manifestation  of  the 

A  A 

discrepancy  between  hQ(k)  and  h(k)  for  small  k.  It  emphasizes  the  convenience 
of  considering  the  spectral  representation  of  the  structure  rather  than  g(r) 
itself.  While  the  attractive  forces  appreciably  contribute  to  only  a  small 
part  of  the  spectrum,  kc^n,  the  attractions  and  repulsions  contribute  (at 
low  and  moderate  densities)  to  g(r)  for  all  r>j.  Thus,  it  is  difficult  to 
say  what  type  of  interaction  produces  a  certain  effect  in  g(r).  Our  results 
indicate  that  this  difficulty  is  not  encountered  when  the  spectrum  is 
considered. 

Finally,  as  a  further  check  on  the  accuracy  of  the  approximation 
gQ(r)  =  g(r)  at  high  densities  we  replace  g(r)  by  g0(r)  in  the  pressure 
(virial)  equation 


Bp/ p  *  1  ♦  (Bo/6)^r(3w(r)/5r)g(r)di 
and  in  the  energy  equation  for  the  excess  internal  energy 
8AE/N  «  (Bp/2)  /w(r)g (r)dr 


(6) 


(7) 


Once  again,  the  calculations  were  performed  within  the  context  of  our 
approximate  treatment  of  the  reference  system.  A  few  of  our  results  are 
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presented  in  Table  I  and  compared  with  those  obtained  from  molecular  dynamics 
and  the  numerical  solution  of  the  Percus-Yevick  equation  for  the  Lennard- 
Jones  system.^  These  results  arc  representative  of  those  found  when  a  far 
more  extensive  comparison  with  molecular  dynamics  is  made  at  several  high 
densities  (pa^2. 65).^ 

It  is  our  opinion  that  we  have  demonstrated  the  validity  of  the 
hypothesis  stated  at  the  beginning  of  this  report.  In  so  doing  we  have  also 
discovered  a  very  simple  and  accurate  theory  for  the  structure  of  dense  liquids 
namely  the  approximation  g(r)  =  gQ(r)  combined  with  Eq.  4.  In  a  future  publi¬ 
cation, ^  our  ideas  will  be  documented  more  thoroughly.  At  that  time  we  will 

also  discuss,  in  detail,  the  relevance  of  our  hypothesis  to  perturbation 

2 

theories  for  liquids. 
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Table  I.  Approximate  and  molecular  dynamics  results  for  the  pressure 

and  internal  energy  as  obtained  from  Eqs.  6  and  7  respectively. 
Columns  3  and  6  give  the  values  found  by  applying  the  approxi¬ 
mation  gQ(r)  =  g(r).  Columns  4  and  7  give  the  molecular 
dynamics  results  (Ref.  8);  the  error  bounds  here  are  ±.05. 

The  results  obtained  by  using  the  solution  of  the  Percus- 
Yevick  equation  for  g(r)  (Ref.  7)  are  given  in  columns  5 
and  8. 


Thermodynamic 

State 

Bp/p 

-BAE/N 

3 

po 

(6c)_1 

This  Kor 

< 

MD 

PY 

This  Work 

MD 

PY 

.85 

.88 

1.82 

1.64 

3.17 

6.77 

6.75 

6.61 

.85 

.786 

1.23 

.99 

2.97 

7.70 

7.70 

7.51 

.85 

.719 

.69 

.36* 

2.82 

8.52 

8.51 

8.28 

.85 

1.128 

2.82 

2.78 

3.57 

5.08 

5.05 

4.98 

*  Molecular  dynamics  results  for  po^=.85,  (Be)"1*. 72,  predict 


8p/ps .43  (Ref.  8) . 


Fig.  1 


Fig.  2 


» 
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FIGURE  CAPTIONS 

A  3 

.  Plots  of  h(k)  for  two  thermodynamic  states:  state  I  is  po  *.5426, 

(8e)_1=.l .326;  state  II  is  po^=.844,  (Be)’*=.723.  The  lines 

represent  Eq.  4';  the  circles  and  triangles  denote  the  molecular 

dynamics  results  (Ref.  1)  for  states  I  and  II  respectively. 

3  -1 

.  Plots  of  g(r)  at  a  high  density:  po  *.85,  (Be)  *.88.  The  line 
represents  the  approximation  gQ(r)  =  g(r);  the  circles  are  molecular 
dynamics  results  (Ref.  1);  the  dashed  line  represents  the  numerical 
solution  of  the  Percus-Yevick  equation  for  the  Lennard-Jones 
potential  (Ref.  7). 
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